Exact and approximate solutions to Schrodinger's equation with decatic potentials 
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Abstract: The one-dimensional Schrodinger's equation is analysed with regard to the existence of 
exact solutions for decatic polynomial potentials. Under certain conditions on the potential's pa- 
rameters, we show that the decatic polynomial potential V{x) = ax^° + bx^ + cx^ + dx''^ + ex^ , a > is 
exactly solvable. By examining the polynomial solutions of certain linear differential equations with 
polynomial coefficients, the necessary and sufficient conditions for corresponding energy-dependent 
polynomial solutions are given in detail. It is also shown that these polynomials satisfy a four-term 
recurrence relation, whose real roots are the exact energy eigenvalues. Further, it is shown that 
these polynomials generate the eigenfunction solutions of the corresponding Schrodinger equation. 
Further analysis for arbitrary values of the potential parameters using the asymptotic iteration 
method is also presented. 
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^ ■ I. INTRODUCTION 

' The problem of finding analytic (polynomial) solutions of linear differential equations has lost some of its interest in 
recent years. The lack of interest is mainly due to the vast development of numerical methods and the fact that much 
of their study has been superseded by more general work based on the theory of Lie algebra [13, HI, [H, H^, [H, [H^, [fil] . 
, ^ , , Recently, the interest of analytic (polynomial) solutions has been renewed by the analysis of confined and un-confined 
quantum systems and also by the search of closed-form solutions to (Schrodin ger- t ype) d ifferential equations with 
^ ; polynomial coefficients P, ill 0, [IMl MM, H, IH-il SI MM,m MM, HI. A reason for such 

^ ■ interest is that in many problems in quantum mechanics, especially those arising from Schrodinger's equation after 
the separation of the asymptotic-behaviour factor of the wave function, there remains a polynomial-type factor in the 
solution lU], il6, JL?, 22, 2^, 26., 21, M, Sl M M, M, M, M, MM, ifiL, 67] . One can further argue that 

' analytic (polynomial) solutions provide a deeper quantitative insight into the physical model under investigation and in 
• , many cases makes the conceptual understanding of physics straightforward and sometimes intuitive 0, [2J, [25l l45l [59| . 
' Moreover, these solutions (if they are available) are valuable tools for checking and improving numerical methods 
introduced for solving complicated physical systems [13] . 



In the present work we study the exact and approximate solutions of the Schrodinger equation with decatic polynomial 
potential [1, i, M, [H, [H 

^. iJV = — T-T + V{x)^/j = Ei:, V{x)^ax^" + bx^ + cx^ + dx'^ + ex'^ (a > 0, - oo < a; < oo). (1) 

The im por tance of this type of even-power potentials follows from its relevance to different models of charmonium 
system J25l [59l and its connection with the analysis of anharmonic and double- well oscillators [1, E[2l[,[23,Ei,[l, 
IstI . l45l l49l [Sll [sgl . l62j . We show under certain relations between the potential parameters a, b, c, d, and e, we 
have exact solutions for the Schrodinger equation ([1]), in the sense that ip and E are completely determined. These 
relations are expressed in terms of energy-dependent polynomials that satisfies a four-term recurrence relation. We 
show further that these polynomials generate the corresponding eigenfunctions. Our approach depends on analysing 
the analytic (polynomial- type) solutions of a second-order differential equation 



E«5,fea;5"My'- Ew'"'h/ = 0, (a^,o + «5,o^0), (2) 



\k=a / \fc=o 
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where 06,^,05^^ and r^^k, for all fc, are constants independent of the variable x. Thereby, we give the necessary and 
sufficient conditions for its polynomial solutions y = X)fc=o '^kx'^- In the second part of the present study, we introduce 
the asymptotic iteration method to obtain accurate approximate solutions for the decatic potential ^ with 
arbitrary values of the potential parameters. In section II, we set up the Schrodinger equation for the decatic-power 
potential ([ij and we show it reduces to a differential equation of a type similar to equation ([2]). In section III, we 
provides both the necessary and sufficient conditions for the existence of polynomial solutions of the (general) linear 
differential equation In section IV, we examine the exact and approximate solutions of decatic potential ^ and 
report some numerical results illustrate the usefulness of our approach where we compare our results with the existing 
literature. 

II. ONE-DIMENSIONAL DECATIC-POWER POTENTIAL 

In this section, we consider the one-dimensional Schrodinger equation 

-i:"{x) + {V{x) ~ E)tP{x) ^0, V{x)=ax^° + bx^ +cx^ + dx^ + ex^, - 00 < a; < 00, (3) 
where a, b, c, d, and e are real constants with a > 0. Assume the solution has the form 

^{x) = xix) ■ e~'^(^) ( lim (j){x) = 00) . (4) 

Substituting this expression into equation ([3)) gives a second-order linear differential equation for x(^) 

X''{x) = 2^\x)x'{x) + - (0'(x))' + V{x) - e)x{x). (5) 

Without loss of generality, we may assume, for the nodeless eigenstate ipo{x), that x(a;) = 1. In which case, equation 
([5]) reduces to Riccati's equation 

(j)"{x) ^ E- V{x) + (^'(x))^ or u'ix) ^E^ V{x) + u^{x), where u(a;) = (ji'^x). (6) 



Definition 1 ffSQ], page J^H) By the symbol y^Pjx) , where P{x) is a polynomial of even-degree, we shall mean the 
polynomial part of the expansion of ^JP{x) in a series of descending integral powers of x. 

Theorem 1 (fSOJ. Theorem 1, page J^H) If in du/dx = Ao + u^, = ^0(2;) is a polynomial of even degree, then no 
polynomial other than 



u = ± 



(7) 



can be a solution of ([5]) . // the degree of is odd, there is no polynomial solution of ([5]) . 
By means of this theorem, we obtain for the solution of equation ([3]) 

where xi^) satisfy the second-order linear differential equation 

64^^ ^ + 8^^ -MxW=0. (9) 
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III. NECESSARY AND SUFFICIENT CONDITIONS FOR POLYNOMIAL SOLUTIONS 



In this section, we study a class of differential equations that generalized equation ([9]) and we establish the necessary 
and sufficient conditions for the existence of polynomial solutions. 

Theorem 2 The necessary condition for the existence of polynomial solutions of the differential equation 

+ (05,02;^ + 05,10;^ + 05,22;^ + 05,30;^ + 05, 4X + 05,5)2/' - (TifiX^ + T4,iX^ + T4,22;^ + r4,3a; + T4,i)y = 0, (10) 



where ago + a-l 7^ and at least one of the coefficients aQj,j = 0, 1, . . . , 6 is different from zero, is 

T4,o = n{n - 1) • 06,0 + n ■ 0,5,0 (n = 0, 1, 2, . . . ) 
while the sufficient condition is given by the solution of the following system of linear equations 

\ 



(11) 
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where 



an 


— n{n 


— 1)06,4 + na^ 


,4 — 7-4,4, 
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(13) 



In particular, for a zero-degree polynomial solution of the differential equation (|10p the sufficient and necessary 
conditions are 



T4,0 = 7-4,1 = 7-4,2 = 7-4,3 = 7-4,4 = 0, 



(14) 



while for a first-degree polynomial solution of equation ()10p . we must have r4.o = 05,0 in addition to vanishing the 
four 2 X 2-determinants 



-T4,4 05,5 

-7-4,3 15,4 — 7-4,4 



-7-4,4 15,5 

-7-4,2 05 3 — T4 3 



0, 



-r4,4 05,5 
-7-4.1 05.2 



7-4.2 



0, 



-T4,4 05,5 

-7-4,0 05,1 — 7-4,1 



0. 



(15) 



For a second-degree polynomial solution, we must have r4,o = 2 • O6,o + 2 • 05,0 in addition to the vanishing of the four 
3 X 3-determinants 



-7-4,4 05,5 

-7^4, 3 05,4 

-7-4,2 05,3 



2a6,6 

74,4 2o6,5 + 2o5,5 

7-4,3 2o6,4 -|- 2o5,4 — r4,4 



-7-4,4 05,5 2o6,6 

-74,3 05,4 — 7-4,4 2o6,5 -f 2o5,5 

-7-4,1 05,2 — 7-4,2 2o6,3 -l- 2o5 3 — r4,3 



-7-4,4 05,5 
-7-4,3 05,4 
-7-4,0 05,1 



7-4,4 
7-4,1 



2a6,6 

2o6,5 + 2o5,5 
206,2 + 205,2 



7"4,2 



-7-4,4 05,5 
-7-4,3 05,4 
^5,0 



7-4,4 
7"4,0 



2a6,6 

2o6,5 -J- 2o5,5 

2o6,i -I- 205,1 



7"4,1 



(16) 
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Similar conditions follows for higher order of polynomial solutions. The first few polynomial solutions are given 
explicitly by 

2/0 (a;) = 1, yi(a;) = 1 + — • 

as, 5 

, s , , -a6,6T-4,3 + (a5,5 + a6,5)7"4,4 , 1 a5,57-4,3 " a5,4T4,4 + t|,4 2 /1 

y2(X) — l-\ 7 r ■ 7 r ■ X + - ■ r TT- • X . (17) 

05,5(05,5 + a6,5j + a6,6(-a5,4 + 7-4,4) ^ (05,5(05,5 + 06,5) + 06,6(^05,4 + ^4,4)) 

This theorem can be proved either by means of the Frobenius series solution of ordinary Differential Equations or 
using the asymptotic iteration method [l^ [E^ ■ 

IV. THE ONE-DIMENSIONAL DECATIC POLYNOMIAL POTENTIALS 

For the existence of polynomial solutions of the differential equation it is clear from the criterion PT|) that the 
parameters a, b, c, and d must be related through the equation 

8{5 + 2n)a^ +8a^d- Aabc + b^ =0, n = 0,l,2,... (18) 

There is an important consequence of this criterion that none of the following decatic polynomials 

V{x):^a-x^°, V{x)=a-x^" + e-x^, V{x) ^ a ■ x^" + c ■ x^ , V{x) ^ a ■ x^° + c ■ x^ + e ■ x^ , (19) 

have analytic solutions of the form ^ with xi^) is a polynomial- type expression. Indeed, we can claim out of this 
criterion that none of these potentials (fT9|) is a quasi-exactly solvable model [sHl and we have to rely on perturbation 
or approximation methods to obtain the necessary information about the discrete spectrum of these potentials as is 
usually used in the literature (e.g. 0, i, [H, [H HI, [H |4l|, [H, [el ) . In the following, we shall focus on studying 
the decatic potentials ([3]) using Theorem [2] to obtain the necessary conditions on the potential parameters for exact 
solutions followed by the application of the asymptotic iteration method to analyse the discrete spectrum for arbitrary 
values of the potential parameters. 

A. Exact solutions 

For the necessary condition for even-degree polynomial solutions of Schrodinger equation ([3|) we have, by means of 
equation ([TT]) that 

8{'in + 5)a^ +8a^d- Aabc + b^ =0, n = 0,l,2,... (20) 
while for the sufficient condition we deduce, by means of equations (|T2| - (fT3|) . that 



where the polynomials P2k{E) satisfies the four-term recurrence relation 

P2n+2{E) = {8a'"^E + {An + 1) {b" - Aac) ) • P2n{E) 

+ 2n ■ {2n - 1) • (32(4n - l)a^/^b - 6"* + 8ab^c - IGa^c^ + 64a^e) • P2n-2{E) 

+ 16384 • • n • (n - 1) • (2n - 1) • (2n - 3) • P2n-iiE), (22) 

initiated with Po{E) = 1, P_„(i?) =0, n> 1. The exact (real) eigenenergies E are the roots of equation P2n+2{E) — 
subject to the following additional condition 

(32{4n+3)a^/%-b'^ + 8ab^c-16a^c^+6Aa^eyP2niE)+4096-a^-n-{2n-l)-P2n-2iE) =0, (n 0, 1, 2, . . . ). (23) 

As an example, for the ground state eigenvalue n = 0, we have for arbitrary values of the potential parameters a, b, 
and c that 

1 4ac - 6^ , 1 40a5/2 + b^ ~ iabc 1 96a^/% - b^ + 8ab^c - IGa^c^ , , 

^0=8 37^' ^ = ~s 2 ' e = 5 24 

8 a'^'-' 8 64 a-' 
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On the other hand, the necessary condition for odd-degree polynomial solutions of Schrodinger equation ([3]) , by means 
of equation (fTTT) . is 

8{4:n + 7)ai +8a^d-4:abc + b^ = 0, (n = 0, 1, 2, . . . ) (25) 
under which the exact wavefunctions are given by 

X2n+3 W = L (2fc + l)!.(2V^)3fe ^ ' " - 0, 1, 2, . . . (26) 
where the polynomials P2k+iiE) are given explicitly using the recurrence relation 

P2n+3iE) - (Sa^/^E + (4n + 3) {b^ ~ 4ac) ) • P2„+i(S) 

+ 2n ■ {2n + 1) • (^32(4n + l)a^^^b - b^ + Sab'^c - IQa^c^ + QAa^e^ ■ P2n-i{E) 

+ 16384- -n • (n- 1) • (2n- 1) • (2n + 1) • P2«-3(^), (n = 0, 1, 2, . . . ) (27) 

initiated with Pi{E) = 1 and P_„ = for n > 1. The exact eigenenergies E are the real roots of P2n+3{E) — 0, 
subject to the additional condition 

(32{4n+5)a^/^b-b^+Sab^c-16a^c^+64a^e^ ■P2n+iiE)+4096-a^ ■n-{2n+l)-P2n-iiE) =0, (n 0, 1, 2, . . . ). (28) 

As an example, for the first excited state n = 1, we have for the potential parameters a, 6, and c that 

3 4ac - 6^ , 1 56a5/2 + ;,3 _ 4^^^ ^ 160a^^^b - b^ + Sab^c - 16a^c^ , , 

^1 ^ a 372— ''^^"o 2 '^ = ~7Ur 3 

8 a^l'^ 8 64 

The first few polynomial solutions, n = 0, 1, ... ,4, are reported explicitly in the appendix. In Table HI we report the 
ground-state eigenvalue for the decatic potentials (|3]) for specific values of the potential parameters a, 6, c, d and e. 
Further, in Table HIl we report the first-excited state for particular values of the potential parameters. These results 
confirm and generalize the early finding by Chaudhuri and Mondal [Tsj , using the improved Hill determinant method, 
for the particular potentials: 

V{x) = ^^x' ~ ^x^~x^^ x^\ {e, = I ^o(^) = exp (^-^ . a:^ + i . a:^ - 1 . x^^ , (30) 

y(x) = ^^x^ - x6 - x« H- x^\ [e^ = \ M^) =x-e^p^~^-x' + l-x^~l- .x^)) . (31) 

The ground-state n ^ and first-excited state n — 1 along with the un-normalized wave functions for the potentials 
(pO)) and (|3T|) are displayed in Figure [1] 




FIG. 1: The ground-state n — and first-excited state n — 1 (un-normalized) wave functions for the decatic potential energy 
V{x) — x^^ + x^ + x^ + dx"^ + ex^ along with the exact eigenvalues as are given by pop and (|3ip . respectively. 
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TABLE I: The ground-state energy Eo for a decatic-power potential V{x) = a • x^^ + b ■ + c ■ + d ■ + e ■ x'^ for different 
values of the parameters. Here, /j. and k are arbitrary positive numbers. 
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TABLE II: The first excited state energy Ei for a decatic-power potential. Here, /x and k are arbitrary positive numbers. 



For higher states, we note, using the resuhs reported in the appendix or using Theorem [5] directly, for the potential 

V{x) = + a;8 + - —x^ 

8 

24(l971 -24^/6423)^ - 8(1971 - 24\/6423)i + 24(1971 + 2476423)^ - 8(1971 + 24\/6423)i - 1551 ^ 
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the exact eigenvalue and wave function are given, in the case n = 2, by 
E2^^ (21 + 4(1971 - 2476423)) 3 + 4(1971 + 2476423))^ 



24 



ij2{x) (1 " + 4(1971 - 24V6423))3 + 4(1971 + 24V6423))3 



3 1 • a;^ ] • exp { -- ■ - - ■ x* — — ■ x'^ 
' ' ^ ' 6 8 16 
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Further, for the potential 

Vix) =x^" -x^ +x^ ^ 

8 

64 1^/3(43 + V42T5) - 64^(23^3 + yT405)(513 + 8^4215)3 + 571(513 + 8V42T5); 



we have the exact solutions 

= i ( 7 + 4- ■ (513 + 8V4215)3 



192(513 + 874215): 
52 



33 
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ip2(x) = I 1 - — f 4 + 4- • (513 + 8\/4215)3 ^ 

^ ^ V 16 V 3i ^ ' 3(513 + 8V42T5)): 



(35) 



(36) 



..^).exp(~l..«-i.x4-A.x^). (37) 
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The plots of these potentials along with their exact solutions for n = 2 are displayed in figure [2] 
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FIG. 2: The decatic potentials (|32ll and (|35[l . respectively, and their exact sohitions 



B. Approximate solutions 



The asymptotic iteration method was introduced Ig*] to investigate the solutions of differential equations of the form 



y" = Ao(x)y' + so(a:)2/, (' = ^) 



(38) 



where Aq = Ao(a;) and so = sq{x) are C°°— differentiable functions. A key feature of this method is to note the 
invariant structure of the right-hand side of ([55)) under further differentiation. Indeed, if we differentiate ([55]) with 
respect to a;, we obtain 



y'" = Ai • y' + si • 2/ 

where Ai = Aq + sq + Aq and s\ = Sq + sqXq. If we find the second derivative of equation psp . we obtain 



(39) 



y^^^ = X2 ■ y' + S2 ■ y 

where A2 — X[ + si + AqAi and S2 — s'l + sqAi. Thus, for {n + 1)*'' and (n + 2)*'' derivative of 
have 



and 



= A„_i • y' + s„_i • y 
= A„ • y' + s„ • 2/ 



respectively, where 

A„ = + Sn-i + AoA„_i and s„ = s^_i + soA„-i. 

From dll]) and (|42]) we have 

A„?;'^"+^' - A„_iy^"+^) = 6„y where (5^ = A„s„_i ~ A„_is„ 



(40) 

n = 1, 2, . . . , we 

(41) 

(42) 
(43) 
(44) 



Clearly, from (l44l) if the solution of (|38p . is a polynomial of degree n, then (5„ = 0. Further, if (5„ = 0, then = 
for all n' > n. In the original paper on the Asymptotic Iteration Method (AIM) 16], the following theorem was 
proved. 



Theorem 3 Given Ag and sq in C°°{a,b), the differential equation i38\) has the general solution 



y{r) = exp — / a{t)dt 



C2 + C1 / exp / (Ao(t) + 2a{T))dT dt 



(45) 
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if for some n > 



= a or equivalently Sn = Sn ■ A„-i ~ s„_i • A„ = 0. 



(46) 



The method can apply directly to the differential equation We have considered here, as an example of the 

effectiveness of the method, the following two decatic potentials (see figure [31) 



and 



V{x) = 0.04a;i° + O.STTx^ + 5.5a;^ - 7.5a;^ + 2x^, 
V{x) = OMx^° + 0.877x8 ^ 5 5^6 _ 7 5^4 _ 3^2^ 



(47) 
(48) 




FIG. 3: Tlie graph of the potentials V{x) = OMx^° + 0.877a;^ + 5.5a;'' - 7.5a;'' + 2x^ and V{x) = OMx^° + 0.877a;* + 5.5a;'' - 
7.5a;'' — 2x^ , respectively. 



For these specific values of the potential parameters the AIM sequences A„ and s„, n = 1, 2, 
with 



can be initiated 



Ao(a;) = 2^x^ + -^x 



b ^3 _ (fcf-foc)^ 
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4ac - 52 
8a3/2 



E 



Aabc 



(49) 



64a3 



8a2 



x\ 



and by means of the iteration sequences ([37]), we compute A„(a;) and s„(a;), for n = 1, 2, 3, . . . . The eigenvalues are then 
the roots of the termination condition (|46p . namely Sn{xo;E) — 0. With several symbolic mathematical programs 
available (Maple; Mathematica, etc.), the computation of the roots of this equation, and thus the eigenvalues, by 
means of the iteration method, provided it is set up correctly, is a straightforward calculation, even for the higher 
iteration steps as shown in Table Hill In principle, the computation of the roots of Sn{xo; E) = should be independent 
of the choice of x = Xg, nevertheless, the right choice of Xq, usually accelerates the rate of convergence to accurate 
eigenvalues within a reasonable number of iterations fl3l . For our computation, we have fixed the initial value of x 
as Xo = 0. 



V. CONCLUSION 



We have discussed exact and approximate solutions of Schrodinger's equations with even-power polynomial potentials. 
Our analysis of exact solutions is based on introducing necessary and sufficient conditions of the polynomial solution 
for a class of differential equations with polynomial coefficients. This analysis can serve as available tools to study 
more complicated physical systems implied by the differential equation. Although we confined our discussion in the 
present work to the decatic potential, our approach can be used to analyse different classes of oscillators with higher 
anharmonicities. For the arbitrary values of the potential parameters, we have applied the asymptotic iteration method 
to solve Schrodinger's equation for which the method proves to be extremely effective for finding the approximate 
solutions. We hope that the study presented here encourages further analysis of (both confined and unconfined) 
polynomial potentials. 
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n 


En{a = 0.04, 6 = 0.877, c = 5.5, d = -7.5, e = 2) 


n 


En{a = 0.04, h = 0.877, c = 5.5, d = -7.5, e = -2) 





1.342 322 779 564 646 216i89 





0.135 429 448 914 739 200i87 


1 


5.136 482 243 202 476 766i96 


1 


2.708 255 201 038 304 023i94 


2 


11.304 205 821 188 349 432203 


2 


8.719 471 783 294 745 896199 


3 


19.577 677 092 617 956 963208 


3 


16.684 057 789 758 348 655204 


4 


29.438 055 052 333 889 749215 


4 


26.234 813 143 437 435 2OI213 


5 


40.723 083 272 298 272 62I222 


5 


37.249 467 272 037 347 758234 


n 


En(a = 0.01, b = 0.1, c = 1.0, d = 3.250, e = 12.5625) 


n 


= 0.01, b = 0.1, c = 1.0, d = 3.050, e = 11.5625) 





3.750 000 000 000 000 000f"^'"=* 





3.611 850 704 712 528 53883 


1 


11.653 048 680 350 115 469io4 


1 


11.250 000 000 000 000 OOOf"^'"^* 


2 


20.358 948 672 177 131 878ii3 


2 


19.713 587 721 031 462 373io9 


3 


29.850 701 419 298 223 478i3o 


3 


28.983 552 585 573 622 185i22 


4 


40.098 623 827 649 000 672i39 


4 


39.026 595 688 356 661 175i43 


5 


51.071 414 767 832 192 856i52 


5 


49.808 260 402 594 522 700i46 



TABLE III: The first six energy eigenvalues (in units of 2m = ft = 1) for tlie one-dimensional decatic potentials V{x) = 
ax^^ + bx^ + cx^ + dx"^ + ex^ for various values of parameters, calculated using AIM with xq = 0. The energies are accurate for 
the number of digits reported. The subscript refer to the number of iteration used by AIM. 
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Appendix 

The first few polynomial solutions of the decatic power potentials using the recurrence relations (|22|) amd (j26[) are: 

• For a zero-degree polynomial solution , n = and X2{x) = 1, 

1. P2{E) := 8a^/^E + - iac = 0, 

2. 96a5/26 - fe'' + 8ab^c - IGa^c^ + Gia^e = 0, 

3. 40a5/2 + - iabc + 8a^d = 0. 

• For a first-degree polynomial solution, n = I and Xsi^) = Cq + CiS, we have cq = 0, and 

1. PsiE) := 8a3/2£; ^ 3(52 _ 4^^) ^ 

2. 160a^/26 -b^ + 8a9c - IGa^c^ + 640^6 = 0, 

3. 56a5/2 _,_ 53 _ 4^5^ 8^2^ ^ 

Thus, if ci ~ 1, 

Xz{x) = X. 

• For a second-degree polynomial solution, n = 2 and Xiix) = cq -I- cix -f C2x'^ , we have ci = and 

1. P^{E) := {Sa^/^E + bib"^ -Aac))P2{E) + 2{ma^'%-b'^ + 8ab'^c-lQa^c^ + Ma^e)PQ{E) = 0, 

2. (224aS/26 - 6^ + Safo^c - IQa^c^ + QAa^e)P2{E) + AQQQa^PaiE) = 0, 

3. 72a5/2 _^ 53 _ 4^5^ 8^2^ ^ 

Further, with ci = 0, C2 = —{^a^E + 6^ — 4ac)co/(16a5), and for co — 1 we have 

Xi{x) ^ Po{E) ~ ——X . 
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• For a third-degree polynomial solution, n = 3 and Xsi^) = cq + cix + C2x'^ + c^x'^, we have cq = C2 = and 

further 

1. P^iE) := {8a^/^E + 5{b^ -4ac))P:i{E)+6{m0a^/%-h^ + 8ab^c-16a^c^ + 64a^e)Pi{E) = 

2. (2880^/25 _ 54 _^ gQ^2^ _ ^Q^2^2 _^ 64a3e)P3(i;) + 122880^ = 

3. 88a5/2 + 53 - 4a6c + 80^^ = 

In addition, for ci = 1, C3 = — P3(i^)/(48a3/2), and the polynomial solution is 

• For a fourth-degree polynomial solution, n = 4 and Xei^) = co + cix -|- C2x'^ + csx^ + 042;^, we have ci = C3 = 
and 

1. PeiE) := [Sa^/^E + 9(62 _ 4^ac))Pi{E) + l2{22Aa^l% - 6^ + gafc^c - IBa^c^ + Ma^e)P2{E)+ 
-M536a5/2(40a5/2 -|- 6^ - 4a6c Sa'^d)Po{E) = 

2. (352a5/26 - 6* + Sab^c - IGa^c^ QAa^e)PA{E) + 24576a5P2(£^) = 

3. 104a5/2 + h^- Aahc + Sa'^d = 

In addition, with cq = 1, C2 = tit; and C4 = — — ^, and the polynomial solution is 

16a'*/"' 1536a'' 

„ P2(£^) 9 Pif-B) , 
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